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Identification of Strongly Correlated Spin Liquid in Herbert- 
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Abstract. - Exotic quantum spin liquid (QSL) is formed with such hypothetic particles as 
fermionic spinons carrying spin 1/2 and no charge. Here we calculate its thermodynamic and re- 
laxation properties. Our calculations unveil the fundamental properties of QSL, forming strongly 
correlated Fermi system located at a fermion condensation quantum phase transition. These are 
in a good agreement with experimental data and allow us to detect the behavior of QSL as that 
observed in heavy fermion metals. We predict that the thermal resistivity of QSL under the ap- 
plication of magnetic fields at fixed temperature demonstrates a very specific behavior. The key 
features of our findings are the presence of spin-charge separation and QSL formed with itinerant 
heavy spinons in herbertsmithite. 
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A search for the materials formed with fermionic spinons 
carrying spin 1/2 and no charge is a challenge for a con- 
densed matter physics pQ. A quantum spin liquid (QSL) 
Can be viewed as an exotic quantum state of matter 
composed of hypothetic particles like chargelcss fermionic 
spinons with spin 1/2 and having a possibility to fluctuate 
as in ordinary fluid. The advances in theoretical studies 
of QSLs date back to Anderson's seminal paper [2] which 
ignited hot discussions in the scientific community. The- 
oretical consensus is that the ground state of QSLs is not 
magnetically ordered, including both gapped and gapless 
spin liquids PQ. However, the experimental investigation 
of QSLs is hindered by a lack of real solids where it can 
occur. In other words, no specific model substances with 
QSL has been found yet, although there are few possible 
candidates. This means that the recognition of a perfect 
material with clear QSL realization makes a major chal- 
lenge in modern condensed matter physics. 

The experimental studies of herbertsmithite 
ZnCu3(OH)6Cl2 single crystal have found no evidences of 
long range magnetic order or spin freezing indicating that 
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ZnCu3(OH)6Cl2 is the promising system to investigate 
QSL [SHE]. Presently herbertsmithite ZnCu 3 (OH) 6 Cl 2 
has been exposed as a S = 1/2 Hcisenberg antiferromag- 
net [3] on a perfect kagome lattice and new experimental 
investigations have revealed its unusual behavior , see 
Rcf. for a recent review. In ZnCu 3 (OH) 6 Cl 2 , the Cu 2 + 
ions with S — 1/2 form the triangular kagome lattice, 
and are separated by nonmagnetic intermediate layers of 
Zn and CI atoms. The planes of the Cu 2+ ions can be 
considered as two-dimensional (2D) layers with negligible 
magnetic interactions along the third dimension. 

A frustration of simple kagome lattice leads to a dis- 
persionless topologically protected branch of the spectrum 
with zero excitation energy known as the flat band [T31TPI] . 
In this case the fermion condensation quantum phase tran- 
sition (FCQPT) [IHTSj can be considered as quantum crit- 
ical point (QCP) of the ZnCu3(OH)6Cl2 QSL composed 
of chargelcss fcrmions with S — 1/2 occupying the cor- 
responding Fermi sphere with the Fermi momentum pp. 
As we are dealing with the 3D compound ZnCu3(OH) 6 Cl2 
rather than with an ideal 2D kagome lattice, we have to 
bear in mind that the real magnetic interactions in the 
substance can shift the QSL from the exact FCQPT point, 
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positioning it in front of or behind the QCP. Therefore, 
the actual location has to be established by analyzing the 
experimental data only. It is believed that the 5 = 1/2 
model on the kagome lattice can be viewed as a gapless 
spin liquid pfl-ll2j , while recent accurate calculations point 
to a fully gapped one, see [T7] and Rcfs. therein. Thus, 
it is of crucial importance to test what kind of QSL is 
formed in herbcrtsmithite and determines its low temper- 
ature thermodynamic and relaxation properties. 





ficients. At the same time, at elevated temperatures a 
non-Fermi liquid (NFL) regime emerges separated by a 
transition (or crossover) region. In other words, the can- 
didate substance should exhibit the LFL, NFL and the 
transition regimes similar to the case of heavy fermion 

(HF) metals and 2D 3 He [IMHQ1QS. This mcans that 
QSL plays the role of HF liquid placed into the insulating 
compound. Spin is carried by an electron which becomes 
localized, while the spin is itinerant. Such a behavior re- 
sembles the spin-charge separation observed recently in ID 
system of interacting non- localized electrons [50]. In our 
case the charges is completely localized while the spins are 
itinerant dwelling on the topologically protected branch. 
As a result, we face a breathtaking picture: New parti- 
cles govern the systems properties at low temperatures 
and these arc different from the particles presented in the 
Hamiltonian describing the system. 
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Fig. 1: Panel A: T-dependence of the magnetic susceptibility 
X at different magnetic fields B [6] shown in the legend. The 
values of Xmax and T max at B = 7 T are also shown. Our 
calculations at B — are depicted by the solid curve x(T) oc 
T~ Q with a — 2/3. Panel B: The heat capacity measured on 
ZnCu3(OH)eCl2 at zero magnetic field [5] is shown by squares. 
Solid curve corresponds to our theoretical approximation based 
on the function C = a\T 3 + d2T 1//3 with fitting parameters ai 
and tX2, see eq. Panel C reports the T-dependence of the 
electronic specific heat C/T of YbRli2Si2 at different magnetic 
fields [21] as shown in the legend. The values of (C/T) max and 
Tmax at B = 8 T are also shown. 

To identify unambiguously the compound having an ex- 
otic QSL located at FCQPT, we must analyze a multitude 
of experimental features rather than a single one. At very 
low temperatures and with magnetic field B applied, we 
expect that the possible QSL substances reveal a Landau 
Fermi-liquid (LFL) behavior of their physical character- 
istics, like the magnetic specific heat and susceptibility, 
spin-lattice relaxation rate and the heat transport coef- 



Fig. 2: Schematic phase diagram of ZnCu3(OH)6Cl2- The 
square at the origin shown by the arrow represents QCP near 
which QSP is located. Vertical and horizontal arrows show 
LFL-NFL and NFL-LFL transitions at fixed B and T respec- 
tively. At fixed T the increase of B drives the system along the 
horizontal arrow from NFL regime to LFL one. On the con- 
trary, at fixed B and raising temperatures the system transits 
along the vertical arrow from LFL regime to NFL one. The in- 
set demonstrates the behavior of the normalized effective mass 
Mff versus normalized temperature Tjv as given by eq. ([5]). It 
is seen that temperatures Tjv ~ 1 signify a transition region be- 
tween the LFL regime with almost constant effective mass and 
NFL one, given by J 1-2 / 3 dependence. It is seen from eq. ([5]l 
that the width of the transition region T w oc T oc B. The tran- 
sition region, where Mjy reaches its maximum at T/T max = 1, 
is shown by the arrows and hatched area both in the main 
panel and in the inset. 

We begin with the analysis of the magnetic susceptibil- 
ity X (T) of ZnCu 3 (OH) 6 Cl 2 shown in the panel A of fig. Q] 
It displays an unusual behavior [BJ: at B > 3 T, x(T) has 
a maximum Xmax(T) at some temperature T max (B). The 
maximum Xmax(T) decreases as magnetic field B grows, 
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while T max {B) shifts to higher T reaching 15 K at B = 14 
T. It is seen from the panel A of fig. [I] that %(T) oc T a 
with a = 2/3. The calculated exponent is in good agree- 
ment with the experimental value a = 2/3 ~ 0.66 [6j. The 
observed behavior of x strongly resembles that in HF met- 
als and associated with their proximity to FCQPT ^BlfTB] . 
The specific heat C, arising from the Cu spin system, at 
the lowest explored temperatures, 106 < T < 400 mK, 
follows a linear temperature dependence, C oc T. As it is 
seen from fig. [TJ panel B, for temperatures of a few Kelvin 
and higher, the specific heat becomes C{T) oc T 3 and is 
dominated by the lattice contribution. At low tempera- 
tures T < 1 K, the strong magnetic field dependence of the 
specific heat C suggests that it is predominately formed by 
the specific heat C mag of QSL, C mag = C — aiT 3 , since the 
lattice contribution is independent of B [3H5] . The above 
behavior of \ ( n g- HI panel A) is a visible parallel to that 
of the HF metal YbRh2Si2 observed in measurements of 
C/T and displayed in the panel C. This coincidence be- 
comes evident if we recollect that for HF liquid x oc C/T 
[T5] . It is seen from the panel C that the electronic specific 
heat of YbRIi2Si2 [H| is also strongly dependent on the 
applied magnetic field and we will see below that both the 
specific heat C mag and that of YbRli2Si2 exhibit the same 
behavior. 
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Fig. 3: The experimental data on measurements of \n [S], 
(C/T)jv = M^f [2T] and our calculations of Afjy at fixed high 
magnetic field (high means that it polarizes completely the 
quasiparticles spins) are shown by points of different shape 
and solid curve respectively. It is clearly seen that the data 
collected on both ZnCu3(OH)6Cl2 and YbRli2Si2 merge into 
the same curve, obeying the scaling behavior. In accordance 
with the phase diagram displayed in the panel A, at growing 
temperatures (y ~ 1) the LFL regime first converts into the 
transition one and then disrupts into the NFL regime. This 
demonstrates that the spin liquid of herbertsmithite is close 
to QCP and behaves like HF liquid of YbRli2Si2 in strong 
magnetic fields |25| . 



To study theoretically the low temperature thermo- 
dynamic, relaxation and scaling properties of herbert- 
smithite, we use the model of homogeneous HF liquid 
|15j . This model permits to avoid complications asso- 
ciated with the crystalline anisotropy of solids. In our 
case, similar to the above HF liquid, QSL is composed of 
chargeless fermions (spinons) with S = 1/2 occupying the 
corresponding two Fermi spheres with the Fermi momen- 
tum pp. The ground state energy E(n) is given by the 
Landau functional depending on the quasiparticle distri- 
bution function n CT (p), where p is the momentum and a is 
the spin index. The effective mass M* is governed by the 
Landau equation [T5"1I!Z2] 
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Here we rewrite the quasiparticle distribution function as 
n a {p,T,B) = n CT (p,T = 0,B = 0) + 5n a {p,T,B). The 
Landau amplitude F is completely defined by the fact that 
the system has to be at QCP of FCQPT [T51I23H25]. see 
[23] for details of solving eq. ([T]). The sole role of Lan- 
dau amplitude is to bring the system to FCQPT point, 
where Fermi surface alters its topology so that the effec- 
tive mass acquires temperature and field dependencies. At 
this point, the term 1 /M* vanishes, eq. (HJ becomes homo- 
geneous and can be solved analytically [T5][23]. At B = 0, 
the effective mass, being strongly T - dependent, demon- 
strates the NFL behavior given by eq. (dJ 



M*(T) ~ a T T- 2/3 . 



(2) 



At finite T, under the application of magnetic field B the 
two Fermi spheres due to the Zceman splitting are dis- 
placed by opposite amounts, the final chemical potential 
fj, remaining the same within corrections of order B 2 . As a 
result, field B drives the system to LFL region, and again 
it follows from eq. dTJ) that 



M*{B) ~ a B B 



-2/3 



(3) 



It is seen from eqs. @ and ^ that effective mass diverges 
at FCQPT. At finite B and T near FCQPT, the solutions 
of eq. (P) M*(B, T) can be well approximated by a simple 
universal interpolating function. This interpolation occurs 
between the LFL regime, given by eq. ([3]) and NFL regime 
given by eq. © [T5"ll2"5] . In the case of strongly correlated 
Fermi systems like HF metals and 2D 3 He the thermo- 
dynamic and relaxation properties are defined by the ef- 
fective mass M*, namely \ « {C/T) oc y/T/T\f oc M*, 
where \/T{T oc \ 2 1S the spin- lattice relaxation rate [15]. 
To study the universal scaling behavior of strongly corre- 
lated Fermi system, it is convenient to introduce the nor- 
malized effective mass and the normalized tempera- 
ture TV dividing the effective mass M* and temperature 
T by their maximal values, Afjjj and Tm respectively. The 
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behavior of = M* l~M* M as a function of y = T/T M 
shown in the inset to fig. [2] is independent of the spe- 
cific features of corresponding strongly correlated Fermi 
system, while both M M and Tm are determined by these 



features 



As a result, we obtain 



XN = (C mag /T) N = (Vl/TxT^ = M* N , 



(4) 



where xn, (C/T) at and (y/l/TiT)^ are the normalized 
values of x, C/T and y/\/T\T, respectively. As a typ- 
ical example, the corresponding maximum values Xmax, 
(C/T) max and T max used to normalize the susceptibility 
X and C/T are shown in fig. [TJ the panels A and C. We 
note that our calculations of Mjy based on eq. (UJ do 
not contain any fitting parameters. The normalized effec- 
tive mass M N = M* /M M as a function of the normalized 
temperature y = T/y = T '/Tm reads 



M* N (y) « co 



ciy 



1 



C2y 



8/3 ' 



(5) 



Here Co = (1 + C2)/(l + ci), c\ and c-i are fitting param- 
eters, approximating the Landau amplitude. Since mag- 
netic field B enters eq. ([TJ only in combination Bfis/T, 
we have T max oc B/ib [15U23] . where \xb is the Bohr mag- 
neton. Thus, for finite magnetic fields variable y becomes 



y = T/T N oc T//j, b B. 



(6) 



Subsequently we use eq. §Q to clarify our calculations 
based on eq. ((TJ). The variables B and T enter eq. (fTJ 
symmetrically, therefore eq. ([5]) determines M* as a func- 
tion of B at fixed T [TS]. It follows directly from eqs. 
© and (TTJ) that 

x{T/vbB) T 2 / 3 ~ X (T, B) T 66 cx 2/ 2 / 3 M^( y ). (7) 

Since the magnetization M(T, B) = J x(T, we obtain 
that 

M(T, B) T~ 1/3 ~ M(T, B) T -0 34 (8) 

depends on the only variable y. Equations ((TJ and ((8j 
confirm the scaling behavior of both xT 66 and MT" 0,34 
experimentally established in Ref. [5] 

Now we construct the schematic phase diagram of 
ZnCu 3 (OH) 6 Cl 2 . This is reported in fig. H At T = 
and B = the system is located near FCQPT without 
tuning. In fact, it located in front of FCQPT since C(T) 
demonstrates the LFL behavior characterized by the lin- 
ear T-dependence at low temperatures @] . Both magnetic 
field B and temperature T play role of the control pa- 
rameters, driving the system from NFL to LFL regions as 
shown by the vertical and horizontal arrows. It follows 
from fig. [3] that in accordance with eq. (|4]) the behavior 
of xn coincides with that of (C/T)jy in YbRh 2 Si2. Figure 
SI the right and left panels correspondingly, demonstrates 
that eqs. ^ and ([3]) are in a good agreement with the 
experimental facts. 
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Fig. 4: Left panel: The temperatures T max (_B) at which the 
maxima of \ ( see fig- UJ the panel A) are located. The solid 
line represents the function T max oc aB, a is a fitting parameter, 
see eq. (©■ Right panel: the maxima Xmax of the functions 
x{T) versus magnetic field B (see fig. [TJ the panel A). The 
solid curve is approximated by ^ lr -..('R s i = r/R -2 / 3 
d is a fitting parameter. 



(B) = dB- 2/3 , seeeq. ©, 
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Fig. 5: Left panel: The specific heat C m a g /T of QSL given 
by eq. (J9j) is extracted from measurements of C(B) on 
ZnCu3(OH)eCl2 at different magnetic fields shown in the leg- 
end Right panel: The temperatures T max (_B) at which the 
maxima of C m ag/T (see the left panel) are located. The solid 
line represents the function T max oc B, see eq. ([SJ). 



According to eq. (£[]), in the case of QSL the behavior 
of the specific heat {C mag /T)N must coincide with that of 
Xn- To separate C mag contribution, we approximate the 
general specific heat C(T) at T > 2 K by the function 

C(T) = ai T 3 + a 2 T 1/3 , 



(9) 



where the first term proportional to a\ is due to the lattice 
(phonon) contribution and the second one is determined 
by the QSL when QSL exhibits the NFL behavior as it 
follows from eq. ©. It is seen from fig. [TJ the panel B, 
that the approximation is valid in the wide temperature 
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range. We note that the value of eti is almost indepen- 
dent of 02, the presence of which allows us to achieve a 
better approximation for C. The obtained heat capacity 
C m ag/T = (C - a 1 T 3 )/T is displayed in the left panel of 
fig. [SJ while the right panel B demonstrates the maxi- 
mum temperature as a function of the magnetic field B. 
It is seen that C mag /T oc M* behaves like \ oc M* shown 
in fig. [TJ the panel A. In figs. H2 and [JJ the normalized 
(Cmag/T) n and xn are depicted. It is seen from both 
these figures that the results obtained on different samples 
and measurements [5J[S] exhibit similar properties. As it 
is seen from figs. [3J [5] and [7J that in accordance with eq. 
(01), (Cmag/T)^ ~ \N displays the same scaling behavior 
as (C/T)m measured on the HF metal YbRli2Si2. There- 
fore, the scaling behavior of the thermodynamic functions 
of hcrbcrtsmithitc is the intrinsic feature of the compound 
and has nothing to do with magnetic impurities. This ob- 
servation rules out a possible supposition that extra Cu 
spins outside the kagomc planes considered as paramag- 
netic weakly interacting impurities could be responsible 
for the divergent behavior of the susceptibility at low tem- 
peratures seen from fig. [TJ the panel A. In that case the 
supposition is to lead to explanations of the observed scal- 
ing behavior of x and C/T in magnetic fields shown in figs. 
[31 [6] and [JJ Obviously, it is impossible. 
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To clarify the scaling behavior in that case, we normalize 
both the function \jT\T and the magnetic field. Namely, 
we normalize (1/TiT) by its value at the inflection point, 
and magnetic field is normalized by Bi n f, Bjv = B/Bi n f. 
It follows from cq. g]) that (1/TiT) N = {M* N ) 2 and we 
expect that different strongly correlated Fermi systems lo- 
cated near FCQPT exhibit the same behavior of the nor- 
malized spin-lattice relaxation rate |15j . It is seen from 
fig. [U that both hcrbertsmithite ZnCu 3 (OH) 6 Cl 2 [26] and 
HF metal YbCu5_ x Au x [27] demonstrate similar behav- 
ior of the normalized spin-lattice relaxation rate. As seen 
from fig. [51 at B < Bi n f (or Bn < 1) the normalized re- 
laxation rate (1/TiT)n depends weakly on the magnetic 
field, while at higher fields (1/TiT)n diminishes in agree- 
ment with eq. ©, (1/TiT) N = {M* N ) 2 oc 5~ 4 / 3 . Thus, 
in accordance with the phase diagram shown in fig. [51 
we conclude that the application of magnetic field B leads 
to crossover from the NFL to LFL behavior and to the 
significant reduction in the relaxation rate. 
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Fig. 7: The normalized susceptibility xn and the normalized 
specific heat (C ma g /T)jv of QSL effective mass versus normal- 
ized temperature Tjv as a function of the magnetic fields shown 
in the legends, xn and (C m ag /T)n are extracted from the data 
of [B] and 5], respectively. 



Fig. 6: The normalized susceptibility xn = x/Xmax = and 
the normalized specific heat (Cmag/T) n = of QSL effec- 
tive mass versus normalized temperature Tjv as a function of 
the magnetic fields shown in the legends, xn is extracted from 
the measurements of the magnetic susceptibility x in magnetic 
fields B [6] shown in the panel A of fig. [T] The normalized 
specific heat is extracted from the data displayed in fig. the 
left panel. Our calculations are depicted by the solid curve 
tracing the scaling behavior of M^. 

Figure [8] displays the normalized spin-lattice relaxation 
rates (1/TiT)n at fixed temperature versus normalized 
magnetic field B^. We notice from fig. [5] that the mag- 
netic field B progressively reduces 1/TiT, and the spin- 
lattice relaxation rate as a function of B possesses an in- 
flection point at some B = Bi n f shown by the arrow. 



As it was mentioned above, QSL plays a role of HF 
liquid framed into the insulating compound. Thus, we 
expect that QSL in herbertsmithite behaves like the elec- 
tronic liquid in HF metals if the charge of an electron 
were zero. In that case, the thermal resistivity w = LT/k, 
w = wo + A W T 2 , of QSL behaves like the electrical re- 
sistivity p — po + A p T 2 of the electronic liquid, since 
A w represents the contribution of spinon-spinon scatter- 
ing to thermal transport, being analogous to the contri- 
bution A p to charge transport by electron-electron scat- 
tering. Here, L is the Lorenz number, k is the thermal 
conductivity of QSL, po and wq arc residual resistivity 
of the electronic liquid and QSL, respectively, and coeffi- 
cient A w - A p oc (M*) 2 [TH[2B]- Thus, we predict that 
the coefficient A w of the thermal resistivity of QSL under 
the application of magnetic fields at fixed temperature be- 



p-5 



V.R. Shaginyan et al. 



haves like the spin-lattice relaxation rate shown in fig. [SJ 
A W (B) oc l/TiT(B) oc (M*{B)) 2 , while in the LFL re- 
gion at fixed magnetic fields the thermal conductivity is a 
linear function of temperature, ft oc T. 
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Fig. 8: The relaxation properties of herbertsmithite versus 
those of HF metals. The normalized spin-lattice relaxation 
rate (1/2iT)jv at fixed temperature as a function of magnetic 
field: Solid squares correspond to data on (1/T\T)n extracted 
from measurements on ZnCu3(OH)eCl2 [26], while the solid 
triangles correspond to those extracted from measurements on 
YbCu5_ x Au x with x = 0.4 [57]. The inflection point where the 
normalization is taken is shown by the arrow. Our calculations 
based on eq. JT} are depicted by the solid curve tracing the 
scaling behavior of (M^) 2 . 



To conclude, we have calculated the thermodynamic 
and relaxation properties of herbertsmithite. Our calcu- 
lations unveil the fundamental properties of QSL, form- 
ing strongly correlated Fermi system located at FCQPT. 
These are in a good agreement with experimental data and 
allow us to detect the behavior of QSL as that observed 
in heavy fermion metals. The key features of our find- 
ings are the presence of spin-charge separation and QSL 
formed with itinerant heavy spinons in herbertsmithite. 
Herbertsmithite represents a fascinating experimental ex- 
ample when the new particles-spinons, non-existing as 
free ones, totally replace the initial particles entering the 
Hamiltonian to dominate the properties at low tempera- 
tures. It is highly desirable to perform the detailed ex- 
perimental studies of the transport properties of herbert- 
smithite. Also, conditions for spinons to form a supcrfluid 
liquid remain to be explored. 
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